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ABSTRACT 


A magnetospheric field model is presented in which the usually 
assumed toroidal ring current is replaced by a circular disk current 
of finite thickness that extends from the tail to geocentric distances 
less than 3 Rg. The drastic departure of this model from the concept 
of the conventional ring current lies in that the current is continuous 
from the tail to the inner magnetosphere. This conceptual change was 
required to account for the recent results of analysis of the 0G0 3 and 
5 magnetic field observations. In the present model the cross-tail 
current flows along circular arcs concentric with the earth and 
completes circuit via surface currents on the magnetopause. Apart from 
these return currents in the tail magnetopause. Mead's (1964) model is 
used for the field from the magnetopause current. The difference scalar 
field, AB, defined as the difference between the scalar field calculated 
from the present model and the magnitude of the dipole field is found 
to be in gross agreement with the observed AB (i.e., the observed scalar 
field minus a scalar reference geomagnetic field). An updated version 
of the AB contours from the 0G0 3 and 5 observations, which is used for 
the comparison, is presented in this paper. Significant differences 
in details exist, however, between the model and the observed results. 
These differences will provide a guide for making modifications in 
the equatorial current system in future models. 



1 . 


Introduction 


Several magnetospheric field models have been proposed and are 
currently in use (e.g., Mead, 1964; Williams and Mead, 1965; Fairfield, 
1968; Olson, 1969). However, recent analyses of the magnetic field 
observations on the 0G0 3 and 5 satellites have shown features that 
are substantially different from the existing magnetospheric models 
(Sugiura et al., 1971; Sugiura, 1972a, 1972b). Two important conclusions 
were drawn from these studies. First, the field depressions in the 
inner magnetosphere must be caused by an equatorial current sheet which 
extends from the tail to geocentric distances of about 3 R^. or less. 
Secondly, there is a +AB region at middle to high latitudes in the 
nightside magnetosphere. The purpose of the present paper is to show 
that gross features of these aspects can be represented by a model 
that includes a magnetopause current, a neutral sheet current, and a 
disk current of which the intensity is uniform from the inner edge of the 
neutral sheet current to 3 Rg and then decreases inward to zero at 2 R . 

Two significant features are incorporated in this model; they are 
(a) return currents of the neutral sheet current via the magnetopause 
surface and (b) the equatorial current sheet which replaces the usual 
toroidal ring current. The point (a) was taken into consideration in 
the model presented by Schield (1969). In the model described below, 
the neutral sheet current flows along circular arcs concentric with 
the earth so that the circular sheet current representing the ring 
current continues from the tail current without any discontinuity. 

The current density in the near-earth tail estimated by Speiser and 
Ness (1967) shows such a circular pattern rather than the usually 
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assutned rectilinear sheet current flowing perpendicular to the earth- 
sun line. 

The model is not a physical one in two respects. First, the 
magnetopause current is represented by the Mead (1964 ) model, and 
hence is self-consistent only with the dipole field and not with the 
combined field of the dipole, the tail current, and the ring current. 
Secondly, the ring current is not based on any particle distribution. 
Regarding the first point, it does not seem advisable to make the 
elaborate calculations required for the determination of the boundary 
current until we have a reasonably good model for other sources. 

Besides, the effect of the magnetopause current is secondary to the 
ring current effect in the context of the present discussion. As to 
the second question, the intent of this paper is to give a first order 
approximation which would provide a guide for future studies. 

2. Magnetospheric Field Model 

The earth's main field . This is represented by a dipole field 
with an equatorial surface field of 0.31 gauss (1 gauss = 10“^ tesla; 
the unit of gamma is used below, Where ly = 10 tesla). The dipole 
axis is perpendicular to the earth-sun line in the present model. 

The magnetopause current . The magnetic field from the magnetopause 
current is represented by the spherical harmonic expansion of the 
scalar potential given by Mead (1964); the coefficients up to m = 5, 
n = 6 in his Table 1 are used. The distance, r^, from the center of 
the earth to the boundary along the earth-sun line is taken to be 10 Rg. 

The tail current . To describe the tail current we use the 
(geocentric) solar geomagnetic coordinate system, i.e. the x-axis toward 
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the sun, the +z-axis along the dipole axis and toward the north dipole 
pole, and the y-axis completing a right-handed system. The tail current 
consists of two parts: (a) an equatorial cross-tail current and (b) 

return currents flowing on the northern and southern magnetopause 
surface (Fig. 1) , The cross-tail current flows parallel to the xy 
plane between z = + z i 30 that its thickness is 2z x . On any plane 
parallel to the xy plane and within this range of z, the current 
density j.(j x , jy) is expressed by a current function Y(x,y) such that 

j x = and j y = -S^/Sx (1) 

Then V.j_ = 0 is automatically satisfied, and lines of equal Y are current 
flow lines. Currents flow along ellipses of the same eccentricity, e, 
if Y satisfies 

(l-e S ) (y/x) 3Y/dx - SY/dy = 0 (2) 

In the model, circular currents are assumed, but the formulas 
are given below for a more general elliptical flow pattern, from which 
circular currents can be obtained by putting e = 0. 

If jy is specified along the x-axis (actually the - x-axis for the 
tail) by 


j y (x,0) = - (dY/dx) x>y „ o * f G0 

then Y along the x-axis is given by 

x 

Y(x,0) = - J*f(x)dx + Yj 

xi 


(3) 


( 4 ) 
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where x x is the position of the earthward edge of the cross-tail 
current, and is a constant. Equation (4) provides the boundary 

condition for (2) . 

Obviously, any function of x 2 + (1-e 2 ) y 2 (=a 2 ) satisfies (2), 
because a is a constant of an ellipse, a being the semi-axis along 
the x-axis. When f(x) is constant, i.e. when the tail current density 
along the x-axis is independent of x, say f (x) = C, then the functional 
form of ^(a 2 ) can be readily determined: 

* - -C [x 2 + (1-e 2 ) y 2 ]^ (5) 

Now, we take C to be a function of z to give a finite thickness 
to the cross-tail current. Then j x , j and the total current density 
j ( = [j x 2 + j y 2 ]^) are given by 

j x - - C (z) (1-e 2 ) y [x 2 + (1-e 2 ) y 2 ]'^ (6) 

j y = C (z) X [x 2 + (1 -e 2 )/ 5 ]"^ (7) 

j = C(z) {[x 2 +(l-e 2 )y* ] / [x 2 + (l-e 2 )y 2 ]}^ (8) 

On the boundary, j (z) is obtained from (8) by putting y = + D (Fig. 1). 

To make the current divergence- free on the boundary surface the magnitude 
of the surface current J g (z) for -z x £ z <: z x is taken to be 

2 

J s^ z) = I | j( z; y = ± D)dz| 

For z < -z x and z > z x , the boundary is a cylindrical surface of 
radius (D 2 + z x 2 ) 2 , and J (z) = const = J (z, ) . 

S g i- 


( 9 ) 
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We assume the following form for C(z): 

C(z) = K cos 2 (rrz/2z 1 ) (10) 

where K is a constant. Thus C(z) = 0 at z = + z 1 , that is, the 
current diminishes to zero at these values of z- z, is taken to be 2 R„ 
in the model. With (10) for C(z), (9) reduces to 

V 2 ’ * <-^ + s 1 " jf") t [x a +(l-e a ) a D 3 ]/[ x 2 + (l- e a )D a ]) % (11) 

The tail current extends from x = -10 to -40 Rg in the present 
model. The magnetic field from the tail current is then obtained by 
the Biot-Savart law through integration over all the currents. This 
also applies to the field from the ring current described below. The 
constant K was chosen such that the field from the tail current is 30y 
at x = -20 Rg, y = 0, and z - 5 Rg, i.e., just above the current slab. 

The ring current . The ring current is represented by a current 
disk having the same thickness, 2z x , as the cross-tail current. The 
outer radius of the disk equals | x x | so that the disk is in contact 
with the cross-tail current slab on the night side. The thickness is 
uniform from the outer edge, taken to be 10 R„, to 3 R™, and inside 

£li -Cj 

the latter distance it decreases linearly, and symmetrically with 
respect to the equator, to become zero at 2 Rg. In the disk, the 
current is circular and concentric with the earth; the current density, 
which is independent of radial distance, varies with z (for -z x ^ z ^ z x ) 
in the same manner as in the tail so that there is no discontinuity in 
the current density at the interface between the disk current and the 
cross-tail current. In expressing the current density, a current 
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function of the same type as for the tail, with e = 0, is used. 

3. Discussions 

It has been demonstrated that the topology of AB is a useful 
concept in the analysis of magnetic field observations in the magneto- 
sphere (Sugiura, 1972a). Here, AB is the observed scalar field minus 
a scalar reference field. The latest version of contour maps of equal 
AB determined from the 0G0 3 and 5 observations for the noon-midnight 
and dawn-dusk meridians for quiet (Kp = 0-1) and slightly disturbed 
(Kp = 2-3) conditions is presented in Figures 2 through 5. Equal AB 
contours provide a simple, but rather sensitive means of quantitatively 
comparing magnetospheric field models with observations. 

Figures 6 and 7 show AB contours, respectively in the noon-midnight 
and dawn-dusk meridians, determined from the model described in the 
preceding section. A comparison of Figures 2 and 6 for the noon- 
midnight meridian indicates a remarkable similarity in general 
character. The agreement between the observed results and the model 
is striking when compared with the AB contours for the Mead-Williams 
model (without a ring current) and those for the Hof fman- Bracken model 
of a ring current (incorporated in the Mead -Williams model), which 
were shown by Sugiura et al. (1971). However, there are significant 
differences between the contours in Figures 2 and 6 in the following 
respects. (i) Magnitudes of the field depressions near the equator 
at close distances (< 4 Rg) are smaller in the model than in the 
observed results. This is especially conspicuous for noon, (ii) +AB 
at middle latitudes on the nightside is greater in the model than the 
observed. (iii) The magnitudes of -AB at high latitudes in the noon 
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meridian are substantially larger in the model than the observed. 

(iv) The +AB values in the outer region of the equatorial magnetosphere 
are slightly larger in the model than those in the observed. 

The nighttime, middle- latitude +AB in Figure 6 is more like that 
in Figure 4 for Kp = 2-3, but then the equatorial -AB values in 
Figure 6 are more different from Figure 4 than from Figure 2. Magnitudes 
of the daytime, high-latitude -AB in Figure 6 are even greater than 
the corresponding values in Figure 4 for slightly disturbed conditions. 
Thus a simple intensification or reduction of the disk current would 
not result in a better agreement between the model and the observations. 

A comparison of Figures 7 and 3 (or 5) for the dawn-dusk meridian 
shows differences that are consistent with those for the noon meridian. 
The main differences lie in: (i) smaller magnitudes of the equatorial 

-AB in the model than the observed, and (ii) the existence of a large 
-AB region over the pole in the model. Though the data coverage is not 
complete, the observed AB over the northern polar cap is probably only 
slightly negative beyond the geocentric distance of about 7 or 8 R„ for 
the Kp = 0-1 condition (Figs. 2 and 3), and is all positive for Kp = 2-3 
(Figs. 4 and 5). 

The differences between the Ab distributions from the model and 
the observations described above appear to imply (a) that the equatorial 
ring (or disk) current in the model is too weak in the inner magneto- 
sphere especially on the dayside, and • (b) that the observed fields at 
high latitudes on the dayside and over the polar cap are not well 
represented by the model. It is noted that the present model for the 
ring current does not include the eastward current due to the pressure 
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gradient at the inner edge of the particle disk, and hence that 
the deficiency in the field depressions in the innermost magnetosphere 
is actually greater than is apparent in the above comparison if the 
effect of the eastward current is taken into account. The existence of 
indentations in the AB contours in Figure 6 in the outer region of the 
dayside equatorial magnetosphere may imply that the particle disk is 
thicker and less concentrated near the equator than in the model in this 
region. Whether or not the magnitudes of -AB at high latitudes on the 
dayside and over the pole can be reduced sufficiently to match the 
observed values by increasing the ring current intensity is not 
obvious. But it seems unlikely that the very large ~Ab (~ -40y) near 
the magnetopause can be reduced substantially by altering the ring 
current intensity at the equator or by increasing r b in Mead's model 
for the magnetopause current. It appears that the theoretical model for 
the magnetopause current does not represent the observed fields at high 
latitudes. The theory does not include the polar cusp, which is an 
important feature observationally (Heikkila and Winningham, 1971; 

Frank, 1971; Fairfield and Ness, 1972). It is likely that the current 
system associated with the polar cusp substantially modifies the overall 
AB distribution at high latitudes including the polar cap magnetosphere. 

Figure 8 gives curves representing the relation between the 
geocentric radial distances of the equatorial crossing of field lines 
and the dipole latitudes of their feet at the earth's surface, determined 
from the present model by field-line tracing in the noon and midnight 
meridians. Corresponding curves for the 0G0 3 and 5 observations are 
also given in Figure 8; see Sugiura et al., 1971 for the details. 
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In spite of the discrepancies between the model and the observations 
described above, the general agreement in Figure 8 is remarkable, 
except for the outer magnetosphere at noon. It is cautioned, however, 
that the 0G0 results (broken lines in Fig. 8) are based on estimates 
from the observations and should still be considered as being provisional. 

Several field lines traced in the noon and midnight meridian planes 
for the present model are shown in Figure 9. The field line leading 
to the neutral point in the noon meridian leaves the earth at a dipole 
latitude between 82° and 82.5°. As an example a field line leaving the 
earth's surface at 70° dipole latitude and at 0600 (or 1800) local time 
is traced in Figure 10; in the lower half of the figure, the projection 
of the field line onto the dipole equator is given, and in the top 
half, the field line is drawn on a curved surface containing it. 

Olson (1972) has recently constructed a magnetospheric field model 
with distributed currents to represent the observed AB distribution. 
According to him, latitudinal cutoffs for low-energy cosmic rays 
determined from his model are in better agreement with recent satellite 
observations of cosmic rays than those calculated from earlier models. 
Details of his model are not available at this writing to enable us to 
make a comparison with the present model. 

4. Conclusions 

The magnetospheric field model presented above has an equatorial 
sheet current which is continuous from the tail to geocentric distances 
less than 3 Rg. This model represents the gross features of the AB 
distribution determined from the OGO 3 and 5 magnetic field observations 
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better than the existing models. The general agreement between the 
model and the observations warrants further improvements of the model 
by replacing the artificial disk current by a current system based on 
an appropriate particle distribution. Since the observational data on 
the particles responsible for the equatorial sheet current are not 
adequate, the discussions presented in the preceding section will be 
useful in selecting particle distributions for future models. 
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FIGURES 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


Figure 8. 


Schematic views of the tail current in geocentric solar 
geomagnetic coordinates. 

AB contours in the noon-midnight meridian plane for Kp = 0-1 
The broken lines at the inner and outer boundaries of the 
contour map indicate border lines beyond which data are 
not available or are presently not considered to be 
reliable due to uncertainties in the orbit accuracies. 

See Sugiura et al, (1971) for the details for the 
derivation of AB. 

AB contours in the dawn-dusk meridian plane for Kp = 0-1. 
(See the caption for Fig. 2), 

AB contours in the noon-midnight meridian plane for Kp = 2-3 
(See the caption for Fig. 2). 

AB contours in the dawn-dusk meridian plane for Kp = 2-3. 
(See the caption for Fig. 2). 

Ab contours in the noon-midnight meridian derived from 
the present model. 

AB contours in the dawn-dusk meridian derived from the 
present model. 

Relationship between the geocentric radial distances of the 
equatorial crossing of field lines and the dipole latitudes 
of their feet on the earth's surface. 
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Figure 9. Field lines for the present model in the noon-midnight 
meridian plane. Field lines are labeled by the dipole 
latitudes at which they leave the earth's surface. 

Figure 10. A field line leaving the earth at dipole latitude 70° 
and at 0600 or 1800 local dipole time: in a curved 

(meridian) surface containing the field line in the upper 

half, and the projection onto the equatorial plane in the 
lower half. 
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